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Abstract 

We investigate Lepton-Flavor Violating (LFV) decays of Higgs to muon-tau in the 
Supersymmetric Economical 3-3-1 (SUSYE331) model. In the presence of flavor mix- 
ing in sleptons {/i, f} and large values ofv/v', the ratio of Br{H — )• T~^fi~) /Br{H — )• 
r^r~) can reach non-negligible values O(10~'^), as in many known SUSY models. 
We predict that for the Standard Model Higgs boson, the LHC may detect its de- 
cay to muon and tauon. We also investigate the asymmetry between left and right 
LFV values of corrections and prove that the LFV effects are dominated by the 
left FLV term, which is 0(10'^) times larger than the right LFV term in the limit 
of small values of \li-p\/TnsijSY ■ The contributions of Higgs-mediated effects to the 
decay r — )• /i/x// are also discussed. 
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1 Introduction 



The experimental evidences of non-zero neutrino masses and mixing [1] have 
shown that the Standard Model (SM) of fundamental particles and interac- 
tions must be extended. Among many extensions of the SM known today, 
the models based on gauge symmetry SU(3)c ® SU(3)i ® U(l)x (called 3-3-1 
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models) [2,3] have interesting features. The model requires that the number 
of fermion families iV be a multiple of the quark color in order to cancel 
anomalies, which suggests an interesting connection between the number of 
flavors and the strong color group. If further one uses the condition of QCD 
asymptotic freedom, which is valid only if the number of families of quarks 
is to be less than five, it follows that is equal to 3. In addition, the third 
quark generation has to be different from the first two, so this leads to the 
possible explanation of why top quark is uncharacteristically heavy (sec, for 
example, [4]). The 3-3-1 models can also provide a solution of electric charge 
quantization observed in the nature [5]. 

In one of the 3-3-1 models [3] three SU(3)l lepton triplets are of the form 
{uiJ, where i/f is related to the right-handed component of the neutrino 
field vi (a model with right-handed neutrinos) . The scalar sector of this model 
requires three Higgs triplets, and it is interesting to note that two Higgs triplets 
has the same U(l)x charge with two neutral components at their top and 
bottom. Giving all neutral Higgs fields a vacuum expectation value (VEV), 
we can remove one Higgs triplet. Hence the Higgs sector of the obtained model 
becomes minimal, and it has been called the economical 3-3-1 model [6]. 

The lepton-flavor is absolutely conserved in the SM. Recently, experiments on 
neutrino oscillations have proved that lepton fiavor is not conserved. It leads to 
motivation on a search for the signals of lepton fiavor violations (LFV) beyond 
the SM. Many versions of the extension of Minimal Supersymmetric Standard 
Model (MSSM) with large tan ^ have been investigated in Higgs LFV decay. 
The interesting here is there exists parameter space that predicts the branching 
ratio of these types of decays are very sizable, enough to be detected by present 
colliders such as CERN Large Hadron Collider (LHC) [7] or International 
Linear Collider (ILC) [8]. For example, the SM [9] predicted that branching 
ratio of if /XT is very suppressed. However, in the beyond SM, this ratio 
can reach large values, more than 10"''. In particular, Refs. [10,11] showed 
that, in the MSSM, BR{H /i+r~) ~ 10"^ if mu/MsusY ~ 10"^ . The 
Minimal Supersymmetric Neutrino Seesaw Models (z/MSSM) [12] predicted 
the branching ratio of heavy Higgs LFV decay is of order 10^^ while that of 
light Higgs LFV decay is of order 10~^. For more discussions in details about 
the Higgs LFV decay, readers are referred to [13] for general LFV framework, 
to [14] for two Higgs doublet models, to [13,14,15,16,17,18,19] for MSSM, and 
i/MSSM and to [20] for httle Higgs models (LTH). 

The MSSM has shown that in the limit of large tan ^5, the radiative corrections 
become non-negligible in many Higgs LFV decay processes. For example, refs. 
[15,17] showed that the ratio Rb/t = BR{H hh)/BR{H tt) can be 
distinguished between the MSSM and non-supersymmetric models. The main 
reason is that the Higgs boson couplings to down-type fermions receive a large 
corrections enhanced by tan^. It leads to many interesting decay processes 
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in quark sector such as b ^ s'~f [21,22]. The large value of tan /3 also leads to 
many interesting effects in the lepton sector, especially when the LFV source 
in sleptons is included. 



Recently, the Supersymmetric Economical 3-3-1 model (SUSYE331) has been 
constructed [23]. Apart from interesting features that mentioned in refs. [23,24,25,26], 
the scalar sector is minimal, and therefore it has been called the economical. 
In a series of works [23,24,25,26], we have developed and proved that the 
non-supersymmetric version [6] and supersymmetric version are consistent, 
realistic and very rich in physics. In the previous work [24], we skip the LFV 
source in the soft sector. However, the model predicts more interesting phe- 
nomenology if there exists LFV source in the soft breaking terms. In this 
paper, we will concentrate on LFV Higgs decays to /xr with the presence of 
misalignment of sleptons {/x, f} and their sneutrinos contained in soft breaking 
terms. In SUSYE331 model, for generating fermion masses as well as canceling 
anomaly, one needs four Higgs triplets. In particular, the "up" p° Higgs gives 
mass for neutrinos and the remain, "down" gives mass for charged leptons 
[23,25] and other Higg s give mass for quarks. The ratio of VEVs, namely ^^lo^ , 
is denoted by tan 7 which is similar to tan/3 in MSMS. Hence, the p° and 
Higgs play very important roles if we consider effects of radiative correction 
in lepton sector in the limit of large tan 7. The corrections may cause many 
non-negligible effects, such as the correction of lepton mass, branching ratio 
of LFV Higgs decay.. .On the other hand, the model 331 is the extension of 
SM based on extended gauge symmetry. Therefore, comparing to MSSM, the 
SUSYE331 model contains new gauge bosons and new Higgses as well as their 
superpartners. Because of appearing of new particles, the number of diagrams 
contributing to LFV Higgs decay in SUSYE331 model is predicted more than 
that in MSSM. It leads to LFV in Higgs decay effected in SUSYE331 may 
be larger than in MSSM. Hence, in this work, we investigate the flavor vio- 
lating Higgs coupling in SUSYE331 model, specially we focus on the {/i, r} 
generations. 



Our work is arranged as follows: In Section 2, we review the particle content 
in SUSYE331 model. The analytic expressions of the Higgs effective couplings 
are studied in Section 3. In Section 4, we study a numerical estimation on 
decay H fir at colliders and compare contribution from the left and right 
LFV radiative corrections into the mentioned decay. In this section, we also 
consider the contributions of Higgs exchange to branching ratio of r ^ 3/^ 
decay. In the last section, we summarize our main results. 
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2 Pcirticle content 



Let us give brief report on the particle content in SUSYE331 model [23]. The 
superfields in the anomaly-free model are given by 

£aL = (Pa,L^?:)^~ (1,3,-1/3), r^i- (1,1,1), a = 1,2, 3 (1) 

QiL^{ui, d,, S')^ -(3,3,1/3), (2) 
ul^, u'l ~ (3*, 1, -2/3), dl^ ~ (3*, 1, 1/3), (3) 

QaL^{da,-Ua,d'a)l^(^,S*,0), a = 2, 3, (4) 

~ (3*, 1, -2/3) , dl^, d!:^ ~ (3*, 1, 1/3) , (5) 
X = (x?, r, xQ^ ~ (1, 3, -1/3), p = pO, - (1, 3, 2/3), (6) 
X' = {x'i,x'^,X2f ~ (1,3*, 1/3), pf = ~ (1,3*, -2/3).(7) 

Here we use some new notations as -ipl = {ipRy = V'lj and exotic quarks are 
denoted by usual quarks with prime-superscripts {u' with the electric charge 
q^i — 2/3 and d' with g^/ = ~l/3). The values in each parenthesis show 
corresponding quantum numbers of the (SU(3)c, SU(3)i:, U(l)x) symmetry. 
In this model, the SU(3)i U(l)x gauge group is broken via two steps: 

SU(3)i ® U{l)x ^ SU(2)^ ® U(l)y U(1)q, (8) 

where the VEVs are defined by 

V2{xf = {u, 0, w) , V2{x'f = {u', 0, ^0 (9) 
V2{pf = (0, V, 0) , V2{pY = (0, v', 0) . (10) 

The vector superfields Vc, V and V' containing the usual gauge bosons are 
given in [23,25]. The supersymmetric model possessing a general Lagrangian 
is studied in [25]. In the following, only terms relevant to our calculations are 
displayed. 
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3 Higgs-muon-tauon effective interactions 



In the SUSYE331 model [23,25], at the tree level, the down-type leptons 
(e, /i, r) only couple to the neutral Higgs {p'^) through the Yukawa interaction 
given by 

>C,// = -^(W,V° + H.c.). (11) 

In general case, Aiaf, 7^ 0, the Lagrangian given in (11) not only provides 
mass for the charged leptons but also gives the source of the lepton flavor 
mixing at the tree level. It means that if the couplings Aia^ 7^ with {a ^ b), 
the LFV processes, such as Higgs^ /xr, must be existed. In this case, our 
theory predicts very large branching ratios of LFV processes which exceed to 
experimental results discussed in [27]. Hence, in the following calculation, we 
skip the Xiab with {a ^ b) in (11) . 

Let us consider another source of LFV which is caused by slepton mixing. 
More details of slepton mixing, one can find in Appendix C. Because of slep- 
ton mixing, the leading effective interactions of leptons with p°, p'° Higgs can 
appcarc at the one-loop order. In this paper, we will concentrate only on the 
couplings of Higgs with {p, r} leptons. 

In order to consider the p, t flavor mixing at the one loop level, first we rewrite 
the original Lagrangian (11) in terms of two component spinor notations which 
are familiar to those in literature, namely 

- Co^^r = {Y^fllflL + YrTlTL) + H.C., (12) 

where = A122/3, = A133/3. 

At the one- loop level, if we skip all of the terms which are proportional to 
except terms contributing to mass of muon, then Yukawa interactions con- 
taining Higgs-lepton-lepton couplings can be divided into two parts: 

• The lepton-fiavor conversing (LFC) part given by 

- AUc = {y,A'; + Y^Af) plpLp'* + Y^A'^Tlnp^* 

+ (y,A'/ + Y^Af) pIplp" + Y^AirlTLp" + H.c, (13) 

• The lepton-fiavor violating (LFV) part given as 

- ACfv = K {^IrlpL + A^rpItl) 

+ n {^irlPL + /.ir^) + H.c, (14) 
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where all of A]f, Af, Ajf', Af', A^, A^', A^, A^', A^ and A^ are the leading 
effective couplings. 

From now on, for convenience, we use notation A to imply any radiative cor- 
rection of couplings appearing in (13) and (14). Note that A is a dimensionless 
function of mass parameters and A^, A|^ are non-zero value even if we assume 
that there is no flavor mixing in slepton sector. We emphasize that A^ is 
one of quantities affecting on many observable quantities such as the ratio of 
branching ratios Br{H — )> bb)/BR{H — )■ rf). The contribution of A(^ to that 
of branching ratios in the SUSY model is studied in [15,17]. The diagrams 
which contribute to all of As are drawn in Appendix A. 

Now let us construct the total effective Lagrangian for Higgs, muon and tauon 
couplings in terms of physical eigenstates. First we write down the whole 
Lagrangian coming from all of Eqs. (12), (13) and (14) in the matrix form 



p^* + H.c.(15) 



where and yi^ are matrices defined by the following formulas: 



yh 



A^p 



A£ 1 + A^' 



yi 



(16) 



with y = YJY, , A-'' = yAlP + Al^ and Af = y + yAY + Af . 



Because of loop corrections, the mass matrix of the /i, r in (15) is no longer di- 
agonal. In order to find the physical eigenstates of muon and tauon, we expand 
the Higgs p and p' around the vacuum expectation values. As a consequence, 
the mixing mass matrix for the muon and tauon are 



- C^ass^Y^v' (plTl]yi\^'^ \+ H.C, (17) 



where 



yi = yi^+t,yi, = {l + A''^ +A% 
^{l + A( + A%)Y,, 
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with 



^'""^ = ^ = W 
1 + A?' + A?i^' 



1 + A?' + A?t^' 
+ A^t, 
1 + A?' + A?i^ 



(19) 



and 



1 



(20) 



It is easy to see that the mixing mass matrix of muon and tauon given in 
(20) is a general matrix. Finding the mass eigenvalues of left-right leptons is 
equivalent to finding a matrix C satisfying: 



c^ylyeC 



vl 
yl 



yl 



(21) 



In our theory, the matrix C can be found in a form 



C 



-sa Ca 



(22) 



where c\ = cos A, s\ = sin A and A is the rotation angle given by 



t2A = tan(2A) 



i + 4-(e^ + 4)' 

In addition, y^ — dmg{y^, yr) in which (|/^, Ur) are defined as follows 



(23) 



2 / 2 2 / I 2 



(24) 



where 



= 4(e^eR + sl)^ + [l + e| - (ej + el) 



(25) 



Note that the mass eigenvalues of muon and tauon are proportional to {y^, y-,-), 
namely 
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(26) 



On the other hand, the mass eigenstates of leptons (/x, r) and (/x'^, r"^) are 
determined from two transformations 







1 

- {Uif 

\ 




L = 


(:) 


= Vi 




y 



ViLl, 



where Ui and Vi have come from (21), namely 



(27) 





■] 












Vt 1 


^Sa 


Ca j 




1 j 



^ Ca sa^ 



(28) 



y -sa Ca J 

Next, we replace 3^;^ in Eq. (15) by a new form deduced from Eq.(18) 
Now we have obtained a new expression of (15) as follows 



- £ = n (1 + A?' + A%) ( /.^ r£ ) X 1^ j 
+ Yr[i^irl)yi, 




(p°*- V'°) + H.c., 



(29) 



In the basis of mass eigenstates of the muon and tauon given in Eq. (27), the 
Lagrangian (29) transforms into 



= n(l + A^' + A%)L'^^ y, Lp"^ 

+ YrL''^{u!y2Vi^)L{p"* - t^p'°) + H.c. 



(30) 



It is needed to emphasize that the first term in Eq. (30) generates only masses 
for muon and tauon while the second creates masses as well as give rise to the 
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lepton flavor mixing. Sources of flavor mixing are two off-diagonal elements of 
the matrix (^7/3^2^/) : 



frrU, T/A 4A^e^ , {cjeL - caSa)Ap 

^ CA5A(Ai6L + A°^6^ - A^6k) - 4(A^eL + 6^A^) ^^^^ 

In the further calculations, we consider of (t^A) ^ 1 but large enough 

(as investigated in MSSM) to cause many interesting effects, and we will com- 
ment more details after some numerical calculations. On the other hand, the 
rotation angle given in Eq. (23) is very small, so we can set ca — 1, sa — 
As a result, Eqs. (24), (31) and (32) can be presented as very simple formulas: 



{ulyi,v^)^^^A''^, (ulyi^v^)^^ ^ ai (33) 

and the above LFV Lagrangian also appears in a simple form: 



- Cpv - n(A^A^V + A^T^A^)(p°* - i,p'°) + H.c. (34) 

Finally, in the mass-eigenstate basis for both lepton and Higgs, we obtain the 
effective LFV Lagrangian: 

- Cfv ^ ^/2y^(A^//V + A^r'^^) (s^s^^s^ae - CaSj^PSase) + H.c, (35) 

where (fiSase 4>Sa36 the Higgs mass eigenstates generated from the mix- 
ing of two original Higgs bosons p° and The expressions of the Higgs mass 
eigenstates were introduced in [25]. They are summarized in the Appendix 
C. The emphasis here is that in the general supersymmetric model there ex- 
ist both the leading interactions of the muon, tauon with neutral scalar and 
pseudo scalar Higgs. However, the SUSYE331 model contains only interactions 
among muon, tauon and scalar Higgs. 

The effective couplings given in (35) are widely investigated for many LFV 
low-energy processes, specially in the MSSM [10,18,28]. In this paper we flrst 
concentrate on some simple aspects of LFV in the SUSYE331 model. In par- 
ticular, we are going to consider the LFV in decays of the scalar Higgs, i.e. 
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— >■ T^/i''^, where $° = (fs^se ^'Sasa- First, we start with studying the 
branching ratios of neutral Higgs decay into muon and tauon. The SUSYE331 
model predicts that the formula of these branching ratios is 

= 2(1 + tan^ 7) (1 Al\^+\ \^) ^ t+t"). 

(36) 

This result is similar to that one given in [10], except the absence of angle of 
mixing among Higgses. In the limit of appropriately large tan 7, the effects 
of LFV in the Higgs decay processes is not to be ignored. Hence, our theo- 
retical prediction is not much different from that of previous results given in 
[10,11,15,16,28]. For details, we will study some numerical calculations for the 
branching ratios indicated by Eq.(36). In our paper, we use the assumption for 
slepton mixing presented in Appendix C. The diagrams giving contributions 
to and A^ are shown in Fig.l. The relevant vertices to our calculation are 
presented in Appendix B. 



Using Feynman rules, we can obtain the expression A^ from the diagrams in 
Fig.l, namely: 



AP_/\P I 



AP 



(37) 



where 



-,/2 



Kb- 



A'' ■ 



Ke 



aP / T r '2 2 ~ 2 \ T t 12 2 ~ 2 \ 



9^ 



247r2 
,2 



flpmxCLSL 
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167r2 

-.2 
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J I 2 2 ~ 2 \ J I 2 2 ~ 2 \ 
T / 2 2-2 \ r/2 2-2 \ 



A'' 



— -A'' 



2-2 - 2 > 
2-2 - 2 
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.12 



^^^^ 2887r2' 
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(a) 



'La 



id) 



(e) I ^0* 



(/) ■ .0* 



1^ P2 P2 



Pi Pi 



{k) 



/ 



:'0 ^ As if 



T p" p 



(0 



1^ 



*0 



'La 

/ 



T 



fx- 



Fig. 1. Diagrams contributing to [(a), (6), (c), (d), (e), (/), (A;)] and [{i), (/) 



Also, the receives contributions from two diagrams (i) and (/) of Fig.l too, 



(39) 



where 



2887r 

I 2 ( J I ri ~ 2 ~ 2 \ J I 12 ~ 2 



T t 12 2 ~ 2 \ T t 12 2 ~ 2 \ 



12 ~ 2 ~ 2 - 



(40) 
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Here we have used some new notations 



vr'^i'l-i 




(41) 



where sl, cl and s/j, cr are deduced from mixing angles for left and right 
handed sleptons (for details, see Appendix C.2). The same relations hold for 
sneutrino sector, with corresponding notations for mixing angles Si,^,Suj^,Cu^ 
and Cujf^. The function l3{x, y, z) is similar to that mentioned in literature [11], 



The analytical results appearing in (38) show that contributions from two 
diagrams (e) and (f) to the always are the same magnitude but opposite in 
sign. Therefore the total contribution of these two diagrams to A£ vanishes. 
On the other hand, results obtained from the Eq.(38) show that if we neglect 
the terms of the slepton mixing, namely sl — sr = Si,^ = Sy^ = 0, the 
amounts collected from the ^ class diagrams given in Fig. 1 are all zero. 
This corresponds to the case of lepton flavor conservation: A£ = A^ = 0. 

We also remind that analytical expressions of other A functions can be 
found in Appendix A. These results demonstrate that the values A^^, A'jf and 
A'^ given by expressions (A.l), (A. 2) and (A. 3) obtained at one loop approx- 
imation, do not vanish even if we assume no mixing of the sleptons. These 
quantities create non-negligible effects of the lepton masses. They are widely 
discussed in many previous papers. Another feature of the SUSYE331 model 
that we would like to remind here: there are two independent sources (Yukawa 
coupling at tree level) to create masses of slepton and neutrino sectors. Hence, 
contributions to LFV corrections come from two independent sources: mixing 
of lepton and sneutrino sectors. We assume that the model contains both LFV 
sources. 

Before coming to numerical computation section, it is necessary to note that 
the formulas of LFV corrections, such as As in this case, have not been estab- 
lished for the SUSYE331 model before. So let us give some general comments 
on the formulas of As which discriminate against those in MSSM versions: 

• At the one loop approximation, the effective couplings As are obtained from 
the diagrams such as those hsted in Figs. (A.l) and (A. 2). We can distin- 
guish two types of diagram which give contribution to As. The first type of 
diagram does not include any Higgsino propagators, for example Fig. 1 (k) 
and (1), and they are known as pure gaugino-mediated diagrams. The sec- 
ond, containing at least one Higgsino propagator like remaining diagrams, is 
Higgsino-mediated type. In general case, each of these kinds of diagram may 



xy \og{x/y) + yz \og{y/z) + zx \og{z/x) 
{x - y){y - z){z - x) 



(42) 
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give main contribution to the As depending on regions of mass parameter 
space. If each Higgsino-mediated diagram gives the dominated contribution 
to As that reach single maximum value. In contrast, each A that gains 
values from pure gaugino-mediated diagrams, A^^ for example, is propor- 
tional to \iXp\. Additionally, we can see the analytic expressions of As given 
in (38), (40) and Appendix A. It is well known in beyond MSSM theories 
[10], all of effective couplings As are obtained from both types of diagrams, 
except A^. In the limit of large values of the dominated contributions 
of As are caused by pure gaugino-mediated diagrams. This conclusion also 
is happened in the SUSYE331. However, in the SUSYE331, there are the 
additional SU (3) l gaugino-mediated diagrams. Hence the values of A^, can 
be changed in comparison with other models. Details of this difference are 
discussed in section 4. 
• The difference between the predictions of the model under consideration 
and other ones due to hypercharge structure of particle content. For ex- 
ample, let us compare our expressions of As with those of As in MSMS 
[10]. All contributions to the As obtained from Fig.l, are proportional to Is 
functions. Rate coefficients in both models are the same level for diagrams 
of Higgsino-mediated type whereas the rate coefficients in the model un- 
der consideration are smaller than that in the MSSM model for diagrams 
of pure gaugino-mediated type. As a consequence, the large contribution 
to the As from the pure gaugino-mediated type will happen if mass pa- 
rameters are large. Furthermore, in this limit of mass parameters, the pure 
gaugino-mediated diagrams are the only source giving contribution to ra- 
diative corrections A^'' of muon mass. It is nature to keep the ratio Yt/Y^, 
at one loop correction, to be the same as those at tree level. This leads to 
the limit of the mass parameters, which does not exceed 10 TeV . 



In the next section we will investigate some numerical results. On that basis, 
we will compare the effects of the LFV origin in the left- and right-slepton 
sectors as well as sneutrino sectors. In order to investigate numerically, we 
are going to use results from [25] such as: g' / g — 3\/2sw ^ _ 0.2312 and 

azi = 128 at the weak scale. 



4 Numerical results 



In this section we firstly discuss some numerical results that relate to any 
signals of LFV decays H — )■ /xr. Let us start with the maximum LFV in both 
left and right sectors, especially si^cj^ = srCr = s^^c^^ = Si,j^c^^ = 0.5. It 
means that we can assign values of mass parameters like: 



13 



5 10 15 20 25 30 2 4 6 8 10 



Fig. 2. |A^p as a function of \fip\/rhji with four different choices of masses ratios: 
1) blue curve-m' = m/j = rhi ; 2) green curve-3m' = rfiR = rfiL ', 3) yehow curve- 
m' = rhji = 2>rhL; 4) red curve-m' = fhpj^ = rhi^j?). Two black lines correspond to 
two values 10"^ and 10"^ of jSOA^p . 





Fig. 3. |A^p as a function of \\Xp\lmL with four different choices of masses ratios: 
1) blue curve~m' = mji = rfiL ; 2) green curve-3m' = rhji = rfiL ; 3) yellow curve- 
m' = ffiL = 3rhji; 4) red curve-m' = itil = m/j/3. A black line corresponds to value 



10-3 of |50A^|2. 



■ 0.2 m? R „ „ 



and 



where ^ are mass parameters used to compare with SUSY mass scale 
^susY- We would like to emphasize that branching ratios of Higgs decays to 
muon and tauon are sizable if tan 7 is large enough. Therefore, in the following 
calculations, we take tan 7 ~ 50. 



The Fig.2 displays the quantity |50A^p as a function of |/ip|/rh,/j while Fig. 3 
displays the jSOA^^p as a function of \^p\/mi where all other relevant parame- 
ters are fixed. Each curve presented in Figs. 2, 3 contains a single peak. All the 
peaks of the curves are obtained at mass parameters at which the contribution 
of the Higgs-mediated diagrams to As are dominated. Corresponding to each 
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curve, there are two regions of mass parameter space separated by deep wells. 
Deep wells, which divide the parameter space into two parts. The first part, 
the mass parameters arc located in the right hand side of deep wells. In this 
region of parameter space, the pure gaugino- mediated type can give main con- 
tribution to As. The second part, the mass parameters are located in the left 
handed side of deep wells at which the dominated contribution to As is ob- 
tained by the Higgs-mediatcd. All of the maximum points of the curves in the 
Fig. 2 and Fig. 3 are reached at \Hp\/ihL^R ~ (9(10~^) and these values depend 
weekly on the changes of values of fhL and rfiR. On the other hand, the max- 
imum values of the A^ is 0(10-^), as concerned in the MSSM [10,12,18,29] 
while the maximum values of the A^ are much smaller than those of A^, 
specifically max(| A^|)^/max(| A^|)^ ~ 10~^. This large difference comes from 
the symmetry of SU (3)^ x U(l)x model. In particular, in the left side of wells 
the main contributions to A^ of SUSYE331 model come from the SU{3)l 
gaugino-mediated diagrams, namely diagrams ((b), (c), (d)) in Fig.l. In con- 
trast, the main contributions to A^ come from only U(l) gaugino-mediated 
diagram. Figs. 2 and 3 also show that both A^ and Ajj are very sensitive 
with the changes of mx, and niR. More details, Fig.4 draws the dependence of 
|A^P/|A^P on lupl/niL, where m' = m\ = rhi = tususy and four different 
fixed values of mR. The maximal and minimal values of the ratio |A^p/| A^p 
on all the curves in Fig.4 have the same value at different values of \iJip\/rhsusY- 
In the parameter region where the Higgs-mediated diagrams give dominated 
contribution to As, the ratio |Ajj.p/|A^p is very small (< 10~^). But in the 
remaining parameters, that ratio is increased. In the limit \lJ^p\/rhsusY ^ 30, 
the ratio |A^p/|A^p reaches a constant value. More general, we can investi- 
gate the influence of ifiR/rfiL on the ratio lAjjp/IA^^p through contour plots 
drawn in Fig.5. On the drawing results showed that |A^p/|A^p < (9(10~^) 
whenever \^p\/msusY ^ 5 and that ratio does not depend too much on the 
ratio friR/rfiL- However, in the limit \iip\/msusY > 7 and rfiR < O.SrfiL, the 
ratio |A^p/|A^p changes very rapidly if small changes and ifiR. It means 
that chirality effects of phenomena relating with A^^ and A^, arc sensitive with 
the change of ratio m^/m^ at large values of /x^. On the other hand, the left 
picture in Fig.5 indicates that when the ratio \^Jip\/TnsusY ^ 7, it will exist 
in some regions of parameter space of ifiR^ifiL at which the contributions of 
left- and right-lepton sectors into the H ^ fir decay process are of the same 
order. In this case, the pure gaugino-mediated diagrams give the dominated 
contribution to both A^ and A^, and also A^^ ( see (A. 2) in Appendix A). Re- 
calling that large values of A^'' can strongly affect directly on the ratio Yp^/Y^-. 
The results presented in Fig.5 again confirm that whenever \Hp\/msusY > 7 
and niR < O.Sm^, the right-lepton sector gives dominated contribution to the 
branching ratio oi H ^ /it decay process. 

Now we investigate more details the region of parameter space where Higgsino- 
mediated diagrams give a dominated contribution. In this region of parame- 
ter space, both A^'' and A^ are much smaller than A^ so we just focus on 
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Fig. 4. |2 as a function of l/ipl/m^ with four different choices of masses ratios: 

1) blue curve-m' = itir = itil ; 2) green curve-3m' = rhu = rriL ; 3) yellow curve- 
m! = rfiL = Srhji; 4) red curve-m' = rhi = rfiji/S. A black line in the left side of 

figure corresponding to the value equals 1. Both black lines in the right side 

lA'' p \ 
of figure presenting are 2 x 10 and 0.1. 




|p,,|/mSUSY 



Fig. 5. Contour plot of j^ji, rfiR/rfiL vs \np\/msusY, where niR = m,y^, 
m\ = rfiL = rhu^ = msusY- The red region corresponds to the values of 
> 0.5. 



m' 



A£. From Fig.6, we can estimate the ratio of Br{H — )■ fiT)/Br{H — t- rr) 
that can reach the order of 10"'^ in the hmit 0.1 < \fJ'p\/MsusY < 6 and 
0.1 < \rhg\/MsusY < 7 where rhg is mass of gauginos. Let us briefly review 
the decay properties of neutral Higgs bosons in the SUSYE331 modeL At the 
tree level, the couplings of neutral Higgs bosons to up-fermions, down-fermion 
are modified with respect to the SM coupling by factors which are given in 
Table 1. 



We assume that all exotic quarks have masses heavier than that of all neutral 
Higgses. It means that the neutral Higgs cannot decay into the exotic quarks. 
The neutral Higgs bosons may decay mainly into the pairs of fermions. This 
prediction depends on the mass of the neutral Higgs. For neutral Higgs (fsase, 
its mass depends on the vacuum expectation values v, v'. So it should be 
predicted SM Higgs with mass smaller than about 130 GeV. Decay of fsase 
to bb and rr are dominated, the branching ratios of 90 percent and 8 per- 
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Fig. 6. Contour plot of BR{H — )• fiT)/BR{H — )• rr), rhg vs Ifipl/msuSY, 
where m' = nix = "nig and rfiR = fhuj^ = friL = rhy^ = msuSY- In 
the left picture, the green and yellow regions correspond to the values of 
BR{H fiT)/BR{H tt) > 0(10-3) . 



Table 1 



Coupling o: 


' neutral Higgs bosons to fermion. 


Particles 


Up-fermion 


Down-fermion 


Exotic up-quark 


Exotic down-quark 


SM Higgs 


1 

Ca 
Sa 


1 

Ca 
Sa 




Sa/ S'j 
Ca 




Ca/ 



cent, respectively. Combined with the results in Fig. 6, the branching ratio 
Br{{psa36 — A*t) is 8 X 10"^ percent. This may be a good signification of new 
physics in the present limits of colliders. For neutral Higgs (psase, it is heavy 
Higgs, the main productions of decay are the the gauge bosons such as W~^W~, 
ZZ,.. .Hence, the branching ratio — )■ /^r is very suppressed. We would like 
to note that the effective interactions of the muon, tauon and Higgs given in 
(35) not only leads to the LFV of Higgs decay process, but also affects the 
other physical processes with lepton-flavor violations. Some of these processes 
which are looked seriously by present experiments are, for instance, r — )■ /i/i/i 
and T — )■ fi'j. Let us apply the effective couplings given in (35) to the r — > yU/i/i 
decay process. In a general way-regardless of the model, the general effective 
Lagrangian describing decay of r — )■ fififi was studied in [30]. However, in 
this work we focus on the effect of the Higgs-mediated LFV interactions on 
the T — )■ fififi decay process. Hence, the four dimensional effective Lagrangian 
which is built through Higgs exchange is formulated by 



C%i^ = -2V2Gpm^m.tan7 + -^^] if^'f^ + m") 

X (A^r> + A^/iV) + H.C.. (43) 

We would hke to remind that the decay process of r — )■ /x/x/i were investigated, 
by [11,30] for examples, in a general model-independent way. The predicted 
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results show that when Higgs exchange effects are much smaller than other 
ones, the ratio i?r(r — )■ 3fi)/BR{T — )■ /i7) becomes constant with a value 
~ C(10~^). Now, we will discuss in more details whether Higgs-mcdiated 
effects can make any significations to the ratio Br{T — )■ ?>ii)/BR{t — )> /i7) in 
the SUSYE331 model. We can divide our results into two cases, namely ^^^^^ 
and </?5^3g Higgs-mediated effects. The results can be written in two respective 
forms: 



1 o mf^ml 



BR(r- ^ /^-/^V-)0j^3g = g tan^ (lA^^ + \^'r\% 



H lo-i 

(44) 

and 



Sa36 8 ' m; 

/tan7\2 /lOOGeV\ 



~ 7 X 10" 



4 



40 / \ m 



/I A/'P I I AP |2\ 

X ( ' ^oi ) (^«) 

These results immediately lead to a consequence: the maximum contribution 
of Higgs exchange processes can be estimated through the formula: 



nr./ + X ^ 11 /tan7\ VlOOGeV\^ 

The values of the branching ratios decrease rapidly corresponding to the en- 
hancement of Higgs masses . We stress that in the model under consideration, 
the (psa36 is identified with the SM Higgs boson and the remain, (pSase^ heavy 
one. Overall, in our model, the SM Higgs-mediated gives larger contribution 
to the branching ratio of the r — )■ fififi decay process than that of heavy Higgs. 
That kind of branching increases if the tan 7 increases. The branching ratio 
estimated in (46) is ~ 10~^^ in the limit of tan 7 ~ 50 and the Higgs mass 
is of the order of 100 GeV. However, the branching ratio of r ^ /x/x/i can be 
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reached at the present hmits of experiment BR{t~ — )■ jjT ^'^jjT) < 3.2 x 10~^ 
[31]. It means that the contribution to r — )■ is suppressed in the hmit 
of tan 7 ^ 50. This result is different from that predicted in the MSSM 
model [16,10,11]. In particular, for the MSSM, the dominant contributions 
to BR{t — > ^J^^J^^J) are induced by the dipole term and the Higgs-mediated 
term at the hmit tan ^ — 50, = 100 GeV. 

Because of sub-dominated contribution of Higgs-mediated to BR{t — )■ yU/x/x), 
the dominated contribution to that branching ratio is still obtained from the 
photon-penguin couplings. This result leads to the values of well-known ratios 
such as 



1^^^^ - 0(10-^) (47) 
BR{t ^ H'j) ^ ^ ^ ^ 

The predicted result is the concerned result given in [16,11,30,32]. 

We emphasize that in the SUSYE331 model, in order to get the dominated 
contribution to the B{t — )■ fJ.fJ.fJ.), the values of tan 7 must be 10^. In this limit 
of tan 7, the result given in (47) is not holden. 



5 Conclusions 



In this paper, we have studied the LFV interactions of Higgs bosons in the 
SUSYE331 model. We have the unique existence of the lepton- number viola- 
tion in the slepton sector at the tree level. On the basis of this assumption 
we have examined the lepton-number violating interactions of Higgs bosons 
at the one-loop level. Specially we have concentrated our study on the LFV 
couplings of Higgs bosons with muon and tauon. The analytical expressions 
of the effective Higgs-muon-tauon couplings are established at the one-loop 
level. One of the features is that the model does not contain the LFV interac- 
tions of neutral pseudo-scalar Higgs bosons. For the neutral Higgs scalars, the 
model contains two types of radiative interactions that violate lepton number, 
namely, (psasefJ'T and (fisasefJ-T- These effective couplings depend on ratios of 
SUSY mass parameters and tan 7. There is an exactly similar to the other 
SUSY models, all LFV couplings are built from two types of diagrams as 
Higgs-mediated diagrams and pure gaugino-mediated ones. Depending on the 
SUSY parameters, each type of diagram gives the main contribution to LFV 
couplings. In this work, we have also studied the branching ratio of the neu- 
tral Higgs decay into muon and tauon. In the limit IfJ.pl/msusY < 7, the ratio 
of BR{H — )■ Tfi)/BR{H — )■ rr) in this region can reach values that can be 
observed by near future experiments and the contributions from both left and 
right LFV sectors to Br{H — > //r) are of the same order. Outside this region 
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the effects of left and right LFV terms mix in different ways in different re- 
gions of mass parameter space. We predicted that for the SM Higgs boson, 
LHC may detect the decay of SM Higgs boson to muon and tauon. For heavy 
Higgs bosons, the branching of LFV decay is very suppressed. We have also 
studied the contribution of Higgs exchange to decay H — > 3//. In the hmit 
tan 7 = 50, the Br{H — >■ 3/x) is very small, out of direct detection of present 
searching of experiment and it leads to predicted results such as the ratio of 
BR{t 3h)/BR{t /X7). 
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A Analytic formulas and diagrams contributing to As 



Let us write down all expressions of As given as following 



I y /27-//2 2~2\i27-//2 2 

""iOStt^^""^ [Ci^3(m ,Hp,mL^) + sMm ,//^, 

9 \ 2 T / 2 2 ~ 2 \ , 27-/ 2 2-2 \' 

„2 



2r/2 2-2 \ i 27-/ 2 2-2 \ 

■^Jsim^, Hp, m^J + s^Jaim^, Hp, m^J 



87r2 



X 



_/2 

y I \ 2 ( 2 T / 12 ~ 2 -2\i 27-/ 12 ~ 2 

^^^^p'rn [cLiciihim ,mL^,mjiJ + Sjih{m ,111^ 

I 2 ( 2j ( 12 ~ 2 ~2\| 27- / 12 ~ 2 ~2 \\' 

+ SL[Cjil3{m ,mL^,mjiJ + Sjil3{m ,mL^,mjiJ) 



+ 



(A.l) 
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'2 ~ 2 ~ 2 \ , r /I2 ~ 2 ~ 2 



^^^pm SlClSrCr 



-him ,mi^^,mjij 



J I 12.-1 ~ Z \ , T t 11- 2. ~ 2 \ 



(A.2) 



9 



9" 

2 



2 T f 2 2 ~ 2 \ , 2 T / 2 2-2 



247r2 

-.2 



9 



167r2 

167r2 



2 T f 2 2~2\i 2r / 2 2~2^ 



27/2 2-2 \| 2r/2 2-2 \ 

4Mml, 11% ml J + clMmlfil, mlj 



2 



87r2 



X 



/2 



(^3(/^p, ^^2) + ^3(/^p, m2^3)) 



P' f3i' '^H2/ ' ^yr^P^ l'L2'' I'm) 
2 /„2 r ^™'2 -2 -2 \ I „2 r ^™'2 -2 -2 



2 2 T f 12-2 - 2 \ X 2 J I 12-2 - 2 ^ 



I 2 ( 2 T / 12-2 -2\i 2 T / 12-2 -2 



(A.3) 



+ c: 



y2 

47r2 
,2 



!^(L-H) 



(/3(/^: 



2 - 2 



(A.4) 



A 2p' ^ g'"^' 

M 1447^2 



X 



h' clCr + h'^SLSR + h! clSr + h' slCr j , ,2 - 2 -2 x 

p. CLCRhim ,mi2,m^J 



, ^UiSii + K^LCr - h' SlCr - H'clSr J f ,2 -2 - 2 X 

+ — ^ SLSRh{m ,m^3,m^3) 

^ T 

-h'p^SLCR + h'.^CLSR - H'^^SlSr + KpCLCR J ( ,2 -2 - 2 



SLCRhim' ^rh^^^rhji^) 



— CLSRh[m ,m^^,mji^) 



(A.5) 
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II I 

A «' A 1 n' 9 



X 



1447r2 



-SLSRh{m ,mL^,mjiJ 



^ h'^SLSR + /i;CLCR - h' SlCr - h' ClSr j , ,2 ~2 ~ 2 N 

H 77 CLCRh[rn .mj^^.mjij 

^ T 

-h' SlCr + h'^CLSR - h' SlSr + H'^lCr w /2 ~ 2 ~2 ^ 

+ — ^r-^^ CLSRh{m ,mL^,mjiJ 

, -h'cLSR + H'^SlCr + h' ClCr - h' SlSr j / ,2 ~ 2 ~2 . 

+ — SLCRhim ,mi^^,mjij 

(A.6) 

In the case of , it also receives contributions from two diagrams (similar to 
diagrams (e) and (f) in Fig.l ) which cancel each other. Therefore we do not 
repeat them in Figs. and A. 3 and A. 4 . Formula of then is 



h'ij_CLCR + h'^SLSR + H'^t-ClSr + h'^^SLCR J ( ,2 ~2 ~2 



X 



Yr 

h'SLSR + h'cLCR - h'sLCR - h'cLSR 



"H-^-L-^u I '"T^i.^u •"^T-^i.^n '"rn^i^-^ti J f 12 ~ 2 ~2 



-h'^SLCR + h'^CLSR - h'^^SLSR + K^ClCr J , ,2 ~2 ~2 



SLCRhim' ^niLs^rfiRj 

T f 12 ~ 2 ~ 2 ^ 



Yr 

+ ■ CLCRh{m ,m^^,m^J 



(A.7) 
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g m 



X 



1447r2 

h'cLCR + h'^SLSR + h' ClSr + H'slCr j , ,2 ~2 ~2 \ 

— TF^ SLCRhim ,mL^,mjiJ 

h'sLSR + h'^CLCR - h' SlCr - h' ClSr -r , ,2 ~2 ~2 \ 

CLSRhim ,mL^,mjiJ 

^ -K^lCr + KclSr - h' SlSr + h' ClCr J I ,2 ~2 ~2 ^ 

+ — Y^ CLCRhim ,m^3,m^J 

-h'cLSR + h'^SLCR + h' ClCr - h'sLSR T , ,2 ~2 ~2 \ 

SLSRhim ,mi^, 771^3) 
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(a) {b) 



(r) (r^) (r) ^ (r^) 



p°* pO* 
(e) .0* (/) 10* 



-I— 



(r) (r<^) (r) (r^) 



/ \ 



Fig. A.l. Diagrams contributing to A^'' ( or 

,pO* 
^--^ --^ Irp 

/ \ 
» / ^ » ^ ^ 

/x Ab 
Fig. A. 2. Diagram contributing to A^''. 



B Lagrangian 



We have denoted by fL,R{fLR) component spinor of the generic mat- 
ter left-handed and right-handed fermion, respectively. The are their 
superpartners which satisfy (/£ = (/i?)^"^, /£ = /^). The four-component 
Dirac spinor can be represented through two-component spinor such as: fi = 
ilJ'L, IJ'r) = ilJ'L, m) and Jl = {fiL, I^r) = {JIl, J^l)- We also emphasize 
that three-left handed leptons contained in the triplet LaL of SU(3)l are 

{falL, fa2L, /aSi)- 
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(a) 



^'0 



(r) (r^ 



'La 



1^ 

ir) 



iv ) 



Fig. A. 3. Diagrams contributing to [(a), (6)], A^/'' [(c)] ( or A^' [(a), (&), (c)]). 



(a) 



A 



B 



Fig. A. 4. Diagrams contributing to A£ [(a)] and A^ [(b)]. 

fall e {UeL, l^f^L, I'tl] = {^^IL, Z^2L, Z^3l}, 
/a2L e {eL,/iL,rL}, 

fa,L e {u:^, U;^, U^^J ^ {U^L, U^,^, Vl^} ^ {{vDl, (z.^^)^, (z/g^)^ (B.l) 

while f^i^ is singlet under SU(3)2,. Conventions for two component spinors 
used in our paper are the same as those given in [11] except the lower index 
L, which is used to distinguish between Dirac spinors / and left-handed Weyl 
spinors fi- 

Next, let us find the interactional vertices relating with our calculation. We 
start to collect related terms from the Lagrangian given in [23,33] into the 
following ones: 

(1) Gaugino and Higgsino mass terms: 



gh ■ 



+ H.C. 



^ 6=1 ^ 

Pp {Ptpr + fp' + PtP2) + H.c. (B.2) 
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where 



(B.3) 



where we have used the results of mass eigenstate states for Higgsinos 
and gauginos given in [24] and [26] . 
(2) Fermion-sfermion-gaugino interaction terms: 



-^{VlXb - L^LXb) + {ITXb - ri^Xe) 



-^{LXLX\-VXLX\), 



(B.4) 



where i = 1,2,. ..,8 is a color index and L = {Ln, L2L, L^l)'^ , L = 
{LiL, L2L, L^lY ■, = {liR, hn, (sr)^ = (^iL) ^2L) ^3lV = {<^R, I^r,tr)'^ = 
(ei, fih and = ([^r, ll^Y = {e%, fl*^, f^, f. In this 

paper we just focus on interactions relating with two fermions, namely 
^aL — {iJ'L,tl}- All interested terms are given as 



r 



fJ'L 



^9' 



.3a/3 " a/2 



a/3 



ig 



>^B + ^={X\-^X\)\i,l 



8a/3 " a/2 



V3 



^(/^i/iiAs - iillitXe) + ^ {iilIilX\ - f^ll^LX\) 



+ (/^lAlA^ - fllfiLX^) 



+ [//^t]. 



(B.5) 



From this, we list the related vertices in Table B.l 
(3) Higgs-Higgsino-gaugino: 



-\{xxXb - ^xXb) + ^(xVAb - x'^x'As) 



2 - - 2 - - 

+ glPpAs - p^pXb) - -{pp'Xb - p'^p'Xb) 

- ^ [MVA^ - ptAVA^ + ^A-^xA^ - X^A-^xA^ 

- ^'A*V'A^ + p'tA-p'A^ - I'A-x'A^ + x'^A-x'A^l (B.6) 
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Table B.l 

Vertices of Icpton-slcpton-gaugino interaction at tree level. 

Vertices of neutral Higgs-Higgsino-gaugino interactions are shown in table 
B.2 

(4) Yukawa interaction terms: 



^llH ■ 



'''' {LaLP%L + LaL~P%j) - ^ {L 



3 



A, 



IIH 



3 L^lILp - '-^e^^\KL)a{L,L)0{p), + H.c. 



(B.7) 
(B.8) 



Our work needs only terms which include leptons p or t, such as: 

Alab 



r — 



3 

Asab 



laLp'Xl + i'^aLpT + Ll^ + KlP2)Kl 



laLptl^bL - laLpl l^bL + ^^uLPi hi - T^aLpikb 



(B.9) 



Prom now we just note that because the conversation of lepton flavor 
in the lepton sector at tree level then Aia^ = with a ^ b and X^cd = 
with c = d. For simplicity, we use new notations: Ye = Am 73,1^^ = 
Ai22/3,y; = A133/3 and Y^^^ = A312/3, = A323/3,K,^ = A313/3. 
Eq.(B.9) now can be written in the common form: 



^iiH — 



PlpIp'^ + {^>^lLp'l + Plp'^ + t>IlP'2 )pI 

TltIp" + {VrLp'r + TLP" + Z^.^P^ )rl 
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Table B.2 

Vertices of the neutral Higgs-Higgsino-gaugino interactions. 
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Table B.3 

Higgsino-lepton-slepton interacions 



- Y^^^llaLpthL - laLpti>bL + ^LptkL - ^aLpthL] 



(B.IO) 



Corresponding vertices are shown in Fig. B.3 
(5) In the soft Lagrangian, the mass term of sleptons is given by 



'/mass 



= ~''^Lab^aL^bL " Rab'' aL^L ~ l^ab^aLP 'bL 



A 



+ habe"^"'{LaL)a{LbL)f:){p)j + -^Hpp*LaLl 

+ ^^^e'^P^{p'*)^{LaL)p{LbL), + H.c. 



c 

bL 



(B.ll) 



here a, 6 are flavor indices {a, 6 = 1,2,3} or a, 6 = {e, //,t} and a, ^,7 
are component indices of SU{2))l. The e°'^'^ is the antisymmetric tensor. 
In this paper we focus on the mixing of slepton Jl and r. This mixing 
makes mass-eigenstate basis of slepton is different from [24]. For more 
detail, please see in Appendix C. The Lagrangian relating with Higgs- 
lepton-slepton interactions has the form 



+ ly,l^pP'*fiL(:tl + Y^^^l^y* {u^^i^U^L - i>^L^rL) + H.C.jB.12) 



Vertices of Higgs- slepton-slepton interactions are hsted in Table B.4. 
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Vertex 


Factor 


Vertex 


Factor 


;;lc jo 
PLTlP 

~c ~ 
P '^f.L'^rL 

p'nri* 




rifLp" 
t^lPlP 
P^v^lKl 

p'pin* 
p"i^;,Ki 


~2^P-Pp 

iYu^rPp 



Table B.4 

Slepton-slepton-Higgs vertices. 



C Mass eigenstates of particles in the SUSYE331 model 



C.l Neutral Higgs 



The physical states of Higgs (mass eigenstates) have been studied in [24] . For 
convenience we review the main results in this appendix. First we expand the 
neutral Higgs components around the VEVs by 

y'^ = ( u+Si+iAi - w+S2+iA2 ] _ ( + t^+gg+iAs + ] 

IT _ ( u'+S-A+iA3 /+ w+Sa+iAa \ n''^ = ( n'- v'+S(,+iAe \ (C, I') 

where {m, w, m', w', t;'}, {S"!, S'a, 5*3, 6*4, S'5, S'g} and {Ai, A2, A3, A4, A5, Ag} 
are VEV, scalar, and pseudo scalar parts of neutral Higgs, respectively. The 
Higgs mass spectrum and the Higgs mass eigenstates given in [24] showed 
that: 

Scalar Higgs: Mass eigenstates of six original scalar Higgs {Si, S2, S3, 5*4, 5*5, Sq} 
are defined as three massless eigenstates {S[^, S'^,ips2A} and three massive 
ones {0S24! V^5a36' '/'Sosfi}- The relations between the original and the mass- 
eigenstate base areLlI: 



^ There are some different definitions for 7 in [23,24,25]. In this paper we use 
notations identifying with those of [23] . 
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^ CHSe 


—S0C0 


— CbCq 


— SaSpSe 


—CaS0S0 




S2 




jj u 




CgSe 









Sr, 














n 


^4 




S0Ce 




-Si3Se 









^5 























[ 








'^a^7 







0524 



where some new notations are defined as follows: 



^ u u „ . „ 

to = tan^ = — = — , ce = cos (J, sg = smtJ, 
w w' 



w 



te = tan ^ = — , cr = cos f3, sq = sin (3, 
w' 



and 



tan 7 



-, s^. = sin7, c-y = cos 7, 



and a is determined through relations: 



tan 2a 



m, 



66a 



m 



-, Cq = COS a, Sa = sma 



33a 



m 



33a 



m; 



18^^ + 9'' 

9g-' + 2g 



/2 



66a 



27 



274 



36a 



9ff' + 2 ^ 
54 



/2 



1 



(^2 + i;'2) (u.'2 + w'^-) 9^2 ^ 2^ 



12 



vw 



54 IcgSyS^I 



(C.2) 



The mass eigenvalues of three physical Higgses 0^24 ) ^Sase (pSase 

+ - V (^33a - ^66a)^ + 4m|6a (C.4) 

+ \/ ("^33a - "^66a)^ + 4m|6a 



m 



33a 



+ m, 



66a 



(C.5) 



Pseudo-scalar Higgs: There are five Goldstone bosons {A^, Aq, A[, A2, ^pa} 
and one massive physical Higgs (pA- Because the (fA does not receive any 
contributions from the A^, Aq and sleptons as well as their sneutrinos only 
couple to p,p', hence there is no coupling of pseudo scalar ipA to muon and 
tauon at the one loop approximation. This is the difference between MSMS 
and our model. 
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C.2 Mass eigenstates of sleptons 



The masses of sleptons in SUSYE331 models were studied in details in [24]. 
In the work, they assumed that lepton numbers are conserved even in the 
slepton sector. This assumption leaded to the absence of mixing terms in 
slepton sector. Our work is interested in studying the source of LFV caused 
by the mixing between slepton fl and f and ignore all other sources of FLV . So 
with two assumptions of R-parity conversation and the small left-right mixing 
in slepton sector, we can base on [24] to write the mass terms of charged 
sleptons in the form: 



At* ^imjljLa + {ml^^jllfL + H.c.) 
+ E ^J*RjRa + [^X,l^*R^R + H.c.) 



(C.6) 



where = {ei, P'L,rL}, Ir, = {ea, /iii, fa} and 



H- 



4' 



mr = Caa — fnt 

iRa 



,'2 



'23 = ^23 



rh%^, = C23 = ml 



6 

_ 1 
~ 6 



^•23' 



■ ^/^02/^03 + ^A22A23(m^ + W^), 

, 9, COS 2/3 ^ 9 cos 27 
' + w^) — 2^ -2^;^ — 

si s^ 



1 



V + «^') (cot2/3-l) - 
9COs2;5 ^ 9 cos 27 



7 

-2^2 



4 L 



*/3 

(cot^ (5 



"7 



(cot^ 7-1) 



1 



2w^ (cot^ 7 - 1 

9 cos 27 , 9 ^ 9, COS 2/3 

si 



^ [-y^ (cot^ 7 - 1) + (xi^ - 2^/;2) (cot^ /3 - l) 
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m:„ in stead of m^^ and m^^^ in the 



Note that here we use notations M^^ and 
soft breaking term of [23]. We assume that the mixing matrix of the ilL,R and 
fj, R slepton masses is given by 



~ 1 ~ 1 

mf, vtlt 



^2 



/ ~2 



(C.8) 



This form is the same as that given in [10,11] ( for detail, see [10], Appendix 
A.l). Hence, the mass eigenstates and eigenvalues of sleptons in our model are 
similar to that in the MSSM [10,11]. In particular, the mass mixing matrix 
of the left handed and right handed of sleptons given in (C.8) produce the 
mass eigenstates such as {Il^Jl-A and {Ir^Jrs}- The corresponding mass 
eigenvalues are {^^iji'^^ig} and {m%.^,m'^,A. 

From now we adopt conventions the flavor states of sleptons are jlL)^L and 
(11, fl while the mass eigenstates are IL2JL3 and Ir^Jr^, respectively. The 
relations between these two kinds of basics are: jiL — clIl2~slIl3, tl = sl^L2 + 
clIl3, with cl = cos6Il, sl = sin6'L; Jll = crIr^ - srIr^, f£ = srIr^ + crIr^, 
with cr — cos9r, sr — sitiOr. The mixing parameters satisfy the following 
relations: 



m 



SlCl 



LjXT 



SrCr 



~ 2 



(C.9) 



C.3 Sneutrinos 



The general Lagrangian which gains masses for sneutrinos is given in [24] as 
follows 




(C.IO) 



where 



(C.ll) 
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and 



( 1 2t^ \ 1 



1 

18' 



2 



Gab = -fl-^^ab I ;^^8 + j + ^ab + ^/^OayUob, 

1 2 

+ Y^A2aA2b(^^^ + + ghcahcbv'^, 



ab 



1 

V2 



(Sab - £ba)t' + -^l^pV'iXsab - Asba) 



(C.12) 



If the LFV happens only in the {z/^, z/^} sector, we can rewrite the non- 
vanishing terms given in (C.12) in more exphcit formulas: 



o2 / I 2t^ \ 1 

= ^aa = \ [H^ + -j=H, - —H,j + Ml + -1,1 

+3^al(.^+z.^)+^.^i:al, 

11 2 

"^Lpr = ^23 = M|3 + -^02/^03 + Y^A22A23(V^ + W^) + -V^A^^2A3c3, 

11 2 

"^Lmt = ^23 = + -/Xo2/^03 + — A22A23(^^^ + u'^) + -^V^ X3c2>^3c3 ■ (C.13) 



Similar to the charged sleptons sector, we denote by {i^^^, i^tl^ ^uri ^tb} the 
flavor eigenstates while by {ul2i i^ls, ^R2i ^rz\ the mass eigenstates. Also, 
notations {^1^.^, 'f^t^.^i "^^^2 ' "^^^as J' denote the mass eigenstates of sneutrinos. 
Here the relations between two bases are: 



Mi 



III 



,2 



III 



I' Rut 



ml — m?. 

^R3 '^R2 



(C.14) 
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